Unitarity as preservation of entropy and entanglement in quantum systems 
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The logical structure of Quantum Mechanics (QM) and its relation to other fundamental principles 
of Nature has been for decades a subject of intensive research. In particular, the question whether 
the dynamical axiom of QM can be derived from other principles has been often considered. In this 
contribution, we show that unitary evolutions arise as a consequences of demanding preservation of 
entropy in the evolution of a single pure quantum system, and preservation of entanglement in the 
evolution of composite quantum systems. 1 



The standard axiomatic formulation of Quantum Me- 
chanics (QM) relies on a set of postulates describing the 
state of a physical system, its time evolution, and the in- 
formation that one can gather about the system by per- 
forming measurements. A specific formulation for joint 
quantum systems - beyond the appropriate extension of 
the postulates - is not regarded in such axiomatic approx- 
imation. However, one of the most genuine and coun- 
terintuitive properties of quantum mechanics - entangle- 
ment - appears only in joint quantum systems. This 
quantum "inseparability" had for many decades ques- 
tioned the logical structure of QM and raised a funda- 
mental philosophical debate JlJ. The most famous ap- 
parent "paradox" aroused by the nonlocal character of 
composite quantum systems, is exposed in the celebrated 
article of Einstein, Podolsky and Rosen Q (the so-called 
EPR paradox) , in which the authors claimed that in any 
plausible physical theory there exist "elements of phys- 
ical reality" and questioned the completeness of QM by 
analyzing entangled states. Afterwards Bell proved in his 
seminal contributions that if one assumes the validity 
of "Einstein locality" , and an underlying hidden variable 
theory, there is an upper limit to the correlation between 
distant events. The claim, in itself, had nothing to do 
with QM, being a general statement about all physical 
theories that have an underlying hidden variable model 
and respects Einstein locality. However, QM predicts 
the existence of certain entangled states whose correla- 
tions (as predicted by QM) violate the upper limit set by 
Bell (Bell inequalities). Violation of the Bell inequalities 
showed that the conjunction of the principle of locality 
(Einstein locality) and the principle of realism (existence 
of a hidden variable model) is incompatible with QM. 
After the discovery of the Bell inequalities, experiments 
were carefully performed to check the predictions of QM 
regarding the correlations in entangled states. However, 
the results of all the experiments performed till date to 
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check for violation of Bell inequalities can be described by 
local hidden variable models, by taking advantage of the 
inefficiencies of experimental apparatuses. These effects 
have come to be known as "loopholes". E. Santos was 
one of the first to point this out Q (see also [j| , and ref- 
erences therein). Curiously enough, in the last decade, 
this puzzling character of entangled states has become 
the seed of an emerging technological revolution, whose 
consequences are hard to overestimate. 
The aim of the present contribution is to establish a link 
between the axioms of QM and the consequences of de- 
manding (i) preservation of entropy in the evolution 
of a single pure quantum system and (ii) preservation 
of entanglement in the evolution of composite quantum 
systems. We shall demonstrate that any of these require- 
ments (either (i) or (ii)), together with other premises, 
retrieves the dynamical postulate of QM, stating that the 
evolution of a quantum system follows the Schrodinger 
equation, i.e. the time evolution must be unitary 6]. In 
both cases, the converse is obviously true. The paper is 
organized as follows. In Sec. we briefly introduce the 
axiomatic formalism of QM and our starting premises 
to derive the dynamical postulate of QM from entropic 
assumptions and probabilistic linearity. We study maps 
(transformations) that preserve disorder (entropy) and 
are linear on mixtures of density matrices. Sec. ITT1 deals 
with preservation of entanglement in composite (bipar- 
tite) pure quantum systems. There, we investigate prop- 
erties of linear maps which preserve entanglement either 
qualitatively or quantitatively. We demonstrate that ev- 
ery linear map that preserves entanglement, has to be 
local, or local after a swap operation. We also show that 
if the preservation of entanglement is quantitative, then 
the map is a local unitary, or a multiple of a local unitary, 
after a swap operation. 



I. SINGLE QUANTUM SYSTEMS 
A. The Axioms of QM revisited 

The well established standard axiomatic formulation 
of QM is based on four different classes of postulates (see 
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e.g. |||) dealing with description, measurements and evo- 
lution of a physical system. The first postulate concerns 
the "static" description of the system, asserting that to 
any physical system there is an associated Hilbert space 
whose dimension corresponds to the number of "degrees 
of freedom" of the physical system. If complete informa- 
tion about the physical system is available, the state of 
the system is represented by an element of this Hilbert 
space (pure state). Otherwise, the state is represented by 
a density matrix, i.e. a self-adjoint, positive semidcfinite, 
normalized matrix acting on the Hilbert space. A second 
class of axioms of QM deals with the measurements and 
their outcomes. There is a third type of axiom dealing 
with indistinguishably of identical particles stating that 
bosons are represented by symmetric states in the corre- 
sponding Hilbert space while fermions by antisymmetric 
ones. These second and third type of axioms will not 
play any role in our considerations below. Finally, there 
exist a fourth type of axiom concerning the dynamics of 
a physical system. It states that the evolution (transfor- 
mation from one state to another) of the system must 
be unitary. Thus, if a system evolves from a state \4>) to 
a state \ip), with \<j>), e H, then it must occur that 
= U where U is a unitary operator, i.e. U'U = I 
(I being the identity operator on the Hilbert space Tt). 
If a physical system evolves from a density matrix p to a 
density matrix p' , then it should happen that p' = U pU* . 

The logical structure of QM, and its relation to other 
fundamental principles of Nature has been for decades 
a s ubj ect of intensive research (see e.g. 0, 0, IS U 
ElllllISEQHEmili)- In particular, the 
question whether the dynamical axiom can be derived 
from other assumptions has been addressed in e.g. [l9j . 
In this paper, we give an alternative derivation to show 
that this dynamical axiom can be derived from requir- 
ing that the evolution preserves disorder (as quantified 
by von Neumann entropy), in the spirit of the second 
law of thermodynamics. Furthermore, we will also con- 
sider a very weak assumption of probabilistic linearity on 
the dynamics of the system. More precisely, we consider 
maps L fulfilling that: 

(i) the disorder of a physical system p is preserved in 
the evolution of the system transforming p to L(p). 
The disorder of a physical system in a state p is 
quantified by its von Neumann entropy S(p), given 
by 

S(p) = -trplog 2/ 9. (1) 

(ii) Probabilistic linearity: If two physical systems, de- 
scribed by states p\ and pi-, evolve to L(p\) and 
L{p 2 ) 1 then a physical system, described by a prob- 
abilistic mixture ppi + (1 — p)p 2 (0 < p < 1) of pi 
and p2, evolves to pL(pi) + (1 — p)L(p2). 

These are the only assumptions that we make on the dy- 
namics of a physical system. We adopt here the "static" 
description of the physical system from Quantum Me- 
chanics. But, by itself, this static description does not 



have any consequences on the dynamics of the system. 
Then we show that the evolution of t he p hysical system 
can be either unitary or anti-unitary 20] . The converse 
of this statement is of course known: Unitary and anti- 
unitary evolutions preserve von Neumann entropy and 
can be defined to be probabilistically linear. As will be 
discussed later, anti-unitary operators cannot describe 
continuous transformations or better to say they cannot 
be continously deformed to unitary operations. Let us re- 
mark here that this result was previously obtained by A. 
Peres [l^ , showing that there is a close relation between 
dynamical evolutions that violate some fundamental ax- 
iom of QM, like unitarity, and those which are forbidden 
by the second law of thermodynamics. We provide here 
an alternative derivation of this result. 



1. The assumption of preservation of disorder 

Let us first briefly discuss the question of quantifi- 
cation of disorder. Under some quite general assump- 
tions, the amount of disorder of a probability distri- 
bution {pi} can be shown to be the Shannon entropy 
H({pi}) = — ^2 i pi\og 2 Pi- This quantifies, on average, 
the amount of information an observer gathers (being 
equal to the disorder that was in the system) when she/he 
gets to know the result of a random variable which is 
distributed according to {pi}- (See Ref. [22] for details.) 
The above observation, along with the fact that the eigen- 
values of a density matrix p are positive and sum up to 
unity (i.e. they can be interpreted as probabilities), led 
von Neumann to define disorder in quantum mechanical 
systems as S(p), given by Eq. QJ, for a physical system in 
a state p. Notice that density matrices can be interpreted 
as statistical mixtures of projectors onto elements of a 
Hilbert space Ji. And S(p) = — Y^i=i Pi 1°S2 Vi-, where 
the pi are the eigenvalues of p, so that the von Neumann 
entropy of p is the Shannon entropy of the eigenvalues of 
p. This is the only reason that we give for our choice of 
von Neumann entropy as the measure of disorder. For an 
axiomatic approach to derive the von Neumann entropy, 
as well as other reasons and advantages of this choice, 
see Ref. [2^. In passing, notice that the introduction of 
the von Neumann entropy does not require that the evo- 
lution is unitary, or anti- unitary |2 II . Note also that the 
assumption of preservation of disorder is not in contra- 
diction with results of statistical mechanics. This point 
is further discussed in Appendix lAl 



2. The assumption of probabilistic linearity 

We now move on to discuss the assumption in item 
(ii). As noted before, density matrices can be interpreted 
as statistical mixtures of projectors onto elements of a 
Hilbert space H. The probabilistic, or statistical inter- 
pretation of density matrices imply that a statistical mix- 
ture of two density matrices is a legitimate density ma- 
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trix, implying by itself a "probabilistic linearity" of the 
dynamics. This probabilistic linearity of density matrices 
really employs the probabilistic (statistical) description 
of states, and that each element of the ensemble (that 
is described by a density matrix) evolves independently. 
This notion of linearity is used e.g. in classical mechanics 
and classical electrodynamics, when we deal with proba- 
bilistic mixtures of states of a system. This is one of the 
reasons, why we say that item (ii) is a "weak" assump- 
tion. Note that the probabilistic linearity is qualitatively 
different than the linearity of the quantum formalism on 
superpositions of vectors on a Hilbert space, which we 
do not assume. Note also that in our case, the evolution 
characterized by the map L transforms a density matrix 
p into a matrix which, in general, is not normalized. 



B. Evolutions that preserve disorder and are 
probabilistically linear are unitary 

The only states that have zero entropy are the pure 
states (one eigenvalue equals unity and the others van- 
ish). Therefore (by item (i)) every pure state has to be 
mapped again to a pure state. Consider two arbitrary dif- 
ferent pure states |0i)(0i| and \<f> 2 )(<f> 2 \- (The bras and 
kets used in this paper are normalized to unity.) Since L 
inevitably transform them to pure states we can write 



L(\<h)(<h\) = <fe|^XiM 



(2) 



Although L does not have to be norm-preserving, L must 
transform density matrices into matrices, which after 
normalization are density matrices. Therefore d\ and d 2 
must be both positive. Otherwise, a statistical mixture of 
the outputs in Eq. J5J , will not be a positive semidefinite 
matrix. 

Consider now the transformation that takes |0i) to 
\ipi) and \<j> 2 ) to \ip 2 ). Unless \<f>i) is orthogonal to \4> 2 ), 
and, since the choice is completely arbitrary, there must 
be a component of \<p 2 ) that is parallel to |</>i), and a 
component that is orthogonal to \4>i). Thus, there is a 
|$) orthogonal to |0i) such that 



h) = Ai|0i) + A 2 |$) 



(3) 



We can draw the phase of A2 into |$), so that we have 
A2 > 0. The same argument applies to the \ip%). In this 
way 

|V*> = Mi 1*), (4) 

with fj>2 > 0, and \^f) orthogonal to Note that 

I Ax ] = \/\ - A| and \p x \ = y/l- juf- The |A»| and |/^| 
are now uniquely defined and depend on the projectors 
|0i)(0i| and Note that A2 and fi 2 cannot vanish. 

Consider now the state 



with < p < 1. The matrix corresponding to p in the 
orthonormal basis {|</>i), 1$)} is 

p+(l-p)|Ai| 2 (l-p)A x A 2 
(l-p)AJAa (l- P )Xj 

where the * denotes a complex conjugation. This matrix 
has the eigenvalues 

si = \{i - yr^AiTvAi), 

s 2 = i(l + ^l-4pA2 + V A 2). 

Due to probabilistic linearity (item (ii)), L maps p onto 
L(p) =pd 1 \i> 1 )(ip 1 \ + (l-p)d 2 \i> 2 )(i> 2 \. (6) 
which has the eigenvalues 

ti = -(a-0), t 2 = -(a + (3), 

where 

a = pdi + d 2 - pd 2 , 



/3 = ^{pdi + d 2 -pd 2 ) 2 + Apdid^Hp - 1). 

For 2x2 matrices with unit trace, their von Neumann 
entropies can be equal if and only if they have the same 
eigenvalues. Since L has to be entropy preserving, the 
ratio between the eigenvalues of p and L(p) must be the 
same, so that 

.,, si h Si t 2 

either — = — , or — = — . 
s 2 t 2 s 2 ii 



A short calculation yields 

\p(di - d 2 ) + d 2 \X 2 



\Jd^d 2 



(7) 



However, \x 2 was introduced in Eq. Q and hence cannot 
depend on the mixing parameter p, which was introduced 
later on. Therefore we must have, for all p, 



S 

Sp 



= 0, 



(8) 



which implies that 



d\ = d 2 



Hence, the two arbitrary pure states \<pi) and \4> 2 ) are 
mapped onto pure states with the same length. Further- 
more (using d\ = d 2 ) it follows that p 2 = \ 2 . But this 
implies that 

Imi| = |Ai|. 

In other words, the modulus of the scalar product is pre- 
served in the evolution: 



p = p\(j> 1 ){<t>i\ + {l-p)\<t> 2 ){(f> 2 \, 



(5) 



|{^l|^)| = |(0l|02)|. 



(9) 
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The conservation of the modulus of the scalar product 
is, however, a very strong condition, since using the 
Wigner's theorem |23 (see Refs. [2|||22| in this regard), 
we obtain now that the transformation that induces 

\<h) - hM, 

\<h) -> |^ 2 ), (10) 

is either linear and unitary or antilinear and anti-unitary. 
Here linearity means that if a transformation induces the 
transitions in Eq. (|10|> . then the same transformation 
induces the transition a\4>\) +b\<f)2) — > a\ipi) +b\ip2)- And 
antilinearity means that the obtained vector is a*\tpi) + 

Continuous transformations cannot be described by 
antilinear operators, as two consecutive antilinear opera- 
tors act as a linear operator. We will disregard the option 
of anti-unitary operators for this reason. The evolution of 
our physical system is not given by Eq. (|10|l . instead, it is 
given by Eq. (0). However, we have shown that d\ = d.2- 
Thus, the evolution on our physical system is unitary, 
up to a constant. As this constant is independent of the 
input state, it is irrelevant. Therefore we have reached 
our goal. The evolutions that respect the second law of 
thermodynamics and are probabilistically linear, are just 
the ones which are postulated in Quantum Mechanics: 
linear (on superpositions of vectors) and unitary. 



respectively m, over the field C. We restrict ourselves to 
linear maps and allow the addition of local ancillas. The 
structure of the maps considered here is the following: 

L:H A B^HA®HA a ®HB®H Ba - (11) 

The above map L could be decomposed into two different 
maps, the first one consisting only in adding local ancillas 
to the initial state, while the second one maps this state 

into the final one: L : H e -» H®HAa®HBa — -* H®HAa® 
TLsa- The first step is obviously entanglement preserving 
since the ancillas are in a product state. Thus for L 
to be entanglement preserving, it is enough to demand 
L' to be entanglement preserving on the image of the 
inclusion. This means that L' acting on any state of the 
form £g) \A a ) <g> Bb) is entanglement preserving, but 
there might be more general states in H®HAa®HBa on 
which L' is not entanglement preserving. 

By dimensional arguments (dim (range L) < 
dimH ab = nm) one can easily realise that there 
exist states in Ha ® Ha<i ® Kb ® Hbo. which do not 
correspond to the image of L, i.e. there exist elements 
of Ha®Hau®Hb ®Hb<i which are not reached by 
mapping the elements of Hab with L. We can, therefore, 
restrict our investigation to a nm-dimensional subspace 
1Z C Ha ® HAa ® Hb ® Hbo, into which all original 
states from H are mapped. Thus, mathematically we 
study linear, entanglement preserving maps of the form: 



II. COMPOSITE QUANTUM SYSTEMS 

In this second part of the paper, we extend our study 
to bipartite (pure) quantum systems. We restrict our 
study to maps that are linear (on superpositions of pure 
states) and preserve entanglement. We will show that 
demanding qualitative preservation of entanglement, i.e. 
mapping separable states onto separable states and en- 
tangled states onto entangled states requires the map to 
be local or the product of a local map and the swap op- 
erator. The stricter requirement of preservation of some 
measure of entanglement requires the local maps to be 
essentially unitaries. 

At this point, it is worth stressing the difference be- 
tween the linearity postulates, assumed in this paper, for 
single and composite quantum systems. In Sec. [I] we 
assumed a probabilistic linearity on the evolution maps, 
which is a rule for mixtures of density operators. In this 
section, on the other hand, evolution maps are assumed 
to be linear on superpositions of pure states. Henceforth, 
the latter version is simply referred to as "linear" maps. 

A. Formalism 

We restrict our discussion to bipartite systems, tra- 
ditionally called Alice and Bob. The Hilbert space as- 
sociated to the bipartite system is denoted by Hab '■= 
Ha®Hb, where Ha/b are Hilbert spaces of dimension n, 



L:H A ®H B ^Tl- (12) 

We will show later that in fact TZ is isomorphic to Hab- 
In the most general case the map L is not necessarily 
norm preserving. To keep track of this fact, as in Sect. I 
we use the "bra" and "ket" notation only for normalized 
states, i.e. = 1. We shall explicitly write down 

any non-unit-length pre-factors (i.e. c\ip)). A map that 
decreases the norm of some states can be physically in- 
terpreted as a device that prepares the final state with 
a non-unit probability. In the remaining cases, the de- 
vice prepares no system for the output. Norm- increasing 
maps are harder to interpret but we allow them for com- 
pleteness. Finally, in our notation, if two states are par- 
allel, i.e. if for a given \ip) and \<j>), there exists a c S C 
such that \ip) = c\(j>), we use the notation \ip) \\ \<j)). Non 
parallel states will be denoted by lf( \cf>). 

B. Entanglement preservation 

To characterize all maps that preserve entanglement, 
we have to first define this property (of entanglement). 
In the following, we distinguish between qualitative en- 
tanglement preservation and the stronger quantitative 
entanglement preservation. Qualitative entanglement 
preservation requires that separable states are mapped 
to separable states and entangled states are mapped to 
entangled states. This type of entanglement preservation 
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leads to the fact that the map must be either local, or 
the product of a local map times the swap operator. 

For a map to be quantitatively entanglement preserv- 
ing, we should demand that there exists an entanglement 
measure E |28|,|29j, such that for all initial states it 
holds E{\i)j)) — E(L(\tp))). It is known that there exists 
exactly one asymptotic measure of entanglement for pure 
normalized bipartite states, which is the von Neumann 
entropy of the reduced system of one of the two parties 
|30|. In the remainder of this paper, we will denote this 
quantity by EQip)). Using naively the entropy function 
also on non-normalized states does not lead to a sensible 
entanglement measure. The entropy function evaluated 
on |c| 2 |a)(a| (the reduced state of the pure product state 
c\ab)) gives — |c| 2 log 2 |c| 2 , which is only for |c| € {0, 1}. 
There are, however, different ways to extend E{\ij))) to- 
wards non-normalized states and in this paper we will 
deal with two possibilities: The first possibility is to nor- 
malize the state, i.e., to use the renormalized measure 
Ei(c\ip)) = E(\ip)). This measure nullifies any physi- 
cal significance that we may try to ascribe to a possi- 
ble norm-change under the map. Alternatively, we may 
keep the norm as a multiplicative pre- factor, i.e., we use 
the probabilistic measure £ l 2(c|i/')) = |c| 2 i?(|V')). In the 
"failed map" interpretation of norm-loss, this measure is 
the average entanglement the device produces in the long 
run. 

Of course, quantitative entanglement preservation im- 
plies qualitative entanglement preservation. But it is the 
quantitative preservation which allow us to demonstrate 
that such a map is, essentially, a local unitary. Before 
going into details, we review for completeness, some well 
established concepts such as the Schmidt decomposition 
for bipartite pure states and the swapping operator. 

Bipartite systems always admit a Schmidt- 
decomposition, i.e. for every pure state \ip) G Hab, there 
exists orthonormal bases {\ai}} of Ha and {\bj)\ of Hb, 
such that = J2l=i A l |a 4 ) ® with A > 31]. The 
Xi are denoted as the Schmidt coefficients and the r as 
the Schmidt rank of the state Hence the Schmidt 
rank is the minimal number of product states needed 
to decompose |?/>), and, by construction, is bounded by 
r < min(n, m). 

Finally, we define here the swap operator as an spe- 
cial map that preserves entanglement. Swapping, in this 
paper, corresponds to exchange or relabeling the two sub- 
systems and is mathematically a map S : Ha ® Hb , — ► 
Ha> ®Hb', where Ha 1 is isomorphic to Hb and Hb> 
is isomorphic to Ha- For a given orthonormal basis 
{|1a), |iu)} of Ha and |m B )} of H B , the 

isomorphisms immediately lead to orthonormal bases 
{\l A ,),...,\m A ,)} of H A - and {|1 B ,), |n B ,)} of H B >- 
We define the swap operator as 

S : |u> ® Lb) ^ \3A') ® |iu'>. (13) 

Of course, a different isomorphism between Ha/b and 
Ha'/b'i that identifies the bases with different 

orthonormal bases leads to a different swap 



operator. But since all local orthonormal bases are con- 
nected by local unitaries, all swap operators are also 
interconnected by local unitaries, that is, up to local 
unitaries, the swap operator is unique. Note here that 
the swap operator considered in this paper is different 
from the usual swap operator considered in the literature, 
which also allows exchanging a subspace of the Hilbert 
space Ha with that of Hb- Such a general swap opera- 
tor can of course increase (or decrease) the entanglement 
between A and B: Simply consider exchanging A" with 
B" (and back) for the state \V~)a'A"\0)b>\0)b" in the 
A' A" : B'B" bipartite cut, where |*") = -i(|01)-|10)). 



C. Qualitative entanglement preservation 

Qualitative entanglement preservation is equivalent to 
the demand that the set of states with Schmidt rank 
r = 1 (product states) and the set of states with Schmidt 
rank r > 1 (entangled states) are both invariant under 
application of the map L. 

Fact 1 A qualitatively entanglement preserving linear 
map has to be of full rank. 

Proof: Suppose that the kernel of a linear qualitatively 
entanglement preserving map L is not empty. Assume 
that there exist a product vector | ai , 6i) which belongs 
to the kernel of L, i.e. L\ ai, b\) — 0. Let us choose a sec- 
ond product vector | 0,2, 62), for which \cl2), | 62) are not 
parallel to \a±), \bi) respectively. Any (nontrivial) com- 
bination of these two product states \ip2), is an entangled 
state of Schmidt rank 2. But unless | 02, 62) is also in the 
kernel, L cannot be entanglement preserving. Thus, if 
there exist one product state in the kernel, then one can 
find a basis of product states all of them belonging to the 
kernel of L, and hence L is the zero map. It is enough to 
show now that there exit no entangled state of Schmidt 
rank 2 in the kernel of L. Assume that there exist one 
entangled state of rank 2 | <j)j) = \\\ e%, f\) + A2 1 e%, /i), 
written in the Schmidt decomposition with Ai > in 
the kernel of L, i.e. L(\ <j>2 )) = 0. We construct 
I </>2 ) = Ai|ei,/i) — A 2 |e 2 ,/2) which is also of rank 2. 
However, (1/2)(| tpfi) + | 0a")) = Ai|ei,/i) and therefore 
L(| 02~)) must be of rank 1, which cannot be since L is 
entanglement preserving map. Therefore, there exist not 
entangled state of rank r = 2 in the Kernel of L. It is 
now obvious to see that if the kernel of L does not con- 
tain any product vector and any Schmidt rank 2 vector, 
then it cannot contain any vector of Schmidt rank r > 2. 
Therefore, we have shown that if L is a non trivial en- 
tanglement preserving map, its kernel is empty and L is 
of full rank. 

We observe that L maps all product states to multi- 
ples of product states. Therefore, L cannot increase the 
Schmidt rank of a pure state: If \<f>) = K\a.ibi}, then 

L\4>) = Yh=i hL\a>ih) = Yh=i A i°i\ a Wi) is a superposi- 
tion of at most r product terms, and the Schmidt rank 
of L\<j)) cannot be larger than r. Furthermore, we know 
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from Fact that L is invertible, that is, for each L that 
is entanglement preserving, there exists an L~ 1 with the 
property that = \<f>). Thus L _1 is also a linear 

qualitatively entanglement preserving map. Thus nei- 
ther L or L^ 1 can increase the Schmidt rank. So, start- 
ing from the weak demand of qualitative entanglement 
preservation we have proved 

Fact 2 The Schmidt rank is invariant under the applica- 
tion of any linear, qualitatively entanglement preserving 
map L. 

We are now ready to state our main result regarding 
qualitative entanglement preservation. 

Result 1 Every linear map that is qualitatively entangle- 
ment preserving has to be either local or a local operator 
times the swap operator. 

The first part of our proof of this result is slightly tech- 
nical, and proceeds by looking at some specific states 
and their image under an entanglement preserving map. 
To do this, we begin by choosing an orthonormal basis 
of Alice's Hilbert space Ha, and an or- 
thonormal basis {|1_b), |ms)} of Bob's Hilbert space 
TLb- Due to the preservation of separability, we know 
that L acts on an element of the product basis |ia, js) 
as: 

IUJ-b) Ci,j\ d id) ® Kj>> ( 14 ) 

where the length Cjj and both tensor factors can a priori 
depend on both input factors \ia) and Since L is 

linear, this evaluation on a basis fixes L completely. Now 
we start by applying L to a state of Schmidt rank 2. Let 
i 7^ k and j ^ I and we look at the transformation 

\1a3b) + lUke) A Ci d \d id ) <g> \e i: j) + c k ,i\d k ,i) ® \e k .i)- 

(15) 

As the right hand is not a product state, we must have 
\dij) K" \dk,i) and \e i: j) jf \e k ,i), i.e., when both indices 
differ, the image vectors cannot be parallel. On the other 
hand, consider the following application of the map to a 
product vector: 

|l A l fl ) + \1 A 2 B ) = \1a) <8> [|1b> + |2 B )] 

c M |d M ) ® |e X) i) + ci )2 |di,2> ® |ei,2) (16) 

It is clear that either (i) \d\ 2) || Mi 1), or (ii) \e± 2) || 
|ei,i). 

It turns out that in case (i), \dij) || Idj^) and \eij) \\ 
while in case (ii), it is the other way around: the 
\d) image vectors only depend on the second index, while 
the \e) image vectors only depend on the first index. 

To prove this, assume that we know that \dij) \\ \di,i) 
for j 7^ / and we look at 

Mb) + Ikytjs) = [\i A ) + |k A )] ® Lis) 

— > CijMij) <g) \e itj ) + c k j\d kd ) <g> |e fcj ), (17) 



for some k ^ i. On the right hand side, the d-factors 
cannot be parallel, since \dij) \\ \di,l) and \dij) jj" \d k ,j) 
as we pointed out before since they differ in both indices. 
Therefore we find that 

31 j= j:\ckj) || \d hl ) =► | eiJ ) || |e fcJ ) Vfc, (18) 

and in an analogous way we can prove that 

3k ? 1 : \e itj ) || \e k>j ) \d itj ) \\ \d hl ) VZ. (19) 

Using [dip) || |di,i) as the starting point, it is now easy 
to show the desired result for case (i): 

\dx, 2 ) II \dx,i) \e it i) || |ei,i) Vz 

=> \di,j) || \di,i) => |ei,j> || (ej-j) 

(20) 

which means that 

lujs) — * Qjlrfij) (8 |ei,j) = Ci,j|di,i)|ejj), (21) 

where c^j may contain extra phase factors as compared 
to Cij. Therefore we have obtained a decomposition of 
L in the form 

L = L A ®L B P 1 (22) 

where L A /b ■ 7~La/b >-> T~Ia/b, with £a|Ia) = |d»,t) and 
is|js) = l e i,j) and P is a global phase/length with P : 
Has !-» Hab with PlUjs) i-> Ci,j|iAjs)- Obviously, the 
map P is diagonal in the chosen product basis, but since 
La ® is is local, it is also clear that P must map all 
product vectors to product vectors in order to fulfill the 
preservation of separability. For case (ii), Eqs. I|18|l and 
(11911 must be replaced by their equivalents with d and e 
swapped. Instead of Eq. (|21|l . we then obtain 

lujs) — * Ci,j\di,j) ® \e iyj ) = Ci,j\dj,j)\e iti ) . (23) 

Since the set {li^js)} spans TLab, and L is of full rank, 
the image of these vectors, i.e. the set {\djj)\eij)}, has 
also to span the whole space 1Z, which is nm dimensional. 
So both sets {\djj)} and {(e^)} have to be linearly inde- 
pendent, and span a m dimensional Hilbert space H.A', 
respectively a n dimensional Hilbert space TLb' , so that 

K = H A > ®H b >- (24) 

In this case, the combined map SoL, where S is the swap 
operator, maps from the space Hab into itself (or more 
specifically in a space isomorphic to TLab), and can be 
decomposed as 

SoL = L A ®L B P, (25) 

where L a \'ia} = (e^), L b \jb) = \dj,j) and P|iAjB) = 
Ci,j IiaJb) • As in the first case, P maps all product vectors 
to product vectors. 

To complete the proof of Result ^ we now only need 
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Fact 3 Every linear map that takes all product states to 
product states, is of full rank, and is diagonal in some 
product basis, has to be local. 

To prove this, we proceed as follows. Let P be of diagonal 
form in the product basis \aibj), that is 



P\aibj) = \ij\aibj). 



(26) 



We know that P maps the product state (J2i a i\ a i)) ® 
Pj\bj)) (where we take the on ^ ^ (3j as arbitrary 
but fixed) to a product state, which can in general be 
written as ^\ . XiDj\aibj). This can only hold if for all i 
and j 



(27) 



Since P is of full rank, all \j ^ 0, and thus, it also has 
to hold that Xii/j ^ and we can write 



Vk Kk fa 



(28) 



The left hand side of eq l|28(l does not depend on i, so 
the fraction Xjk ■— does neither. This gives us a 
separation of the form 

Ay = AjoXjo =: frVj- (29) 

Therefore P — Pa® Pb as stated above, where Pa|o») = 
IJLi\ai) and Pb|6j-) = Vj\bj). 



D. Quantitative entanglement preservation 

In the previous section, we only required a qualitative 
preservation of entanglement. In this section, we will dis- 
cuss the quantitative preservation of entanglement. As 
stated in Sec. Ill Bl a quantitative preservation of entan- 
glement depends on the extension of the von Neumann 
entropy of the reduced system that is used. However, for 
both of the two previously defined entanglement mea- 
sures, E\ and E 2 , we have that 

Result 2 Every local and linear map L that is quantita- 
tively entanglement preserving, has to be a multiple of a 
local unitary. 

We have proved in the previous section that a qualita- 
tively entanglement preserving map is a local map or a 
local map after the application of the swap operator. Ev- 
ery quantitatively entanglement preserving map is also a 
qualitatively entanglement preserving one, so we know 
that L = La <8> Lb ■ Due to the singular value decompo- 
sition of La and Lg we have 



L A = VaDaUa, 



L b = VbDbUb, 



(30) 



where U A , Ub, Va, Vb are unitary matrices, Da is a diag- 
onal matrix in the basis {|1a) 5 l n A)} with the diagonal 



values Ai > A 2 > ... > \ n > 0, and Db is a diagonal ma- 
trix in the basis {Is), |ms)} with the diagonal values 
A*i > A*2 > — > Mm > 0. We define four further sets of 
orthonormal bases by 

k-> :=U^\i A ), \bt) :=UB l \i B ), (31a) 

\ei):=V A \i A ), \fi):=V B \i B ). (31b) 

For the respective proofs we will deal with a special set 
of states defined as 



\i>{c)) = c\a x bi) + |c| 2 K& m », (32) 

where c € [0, 1]. These states transform under L into 



L\ip(c)} = cAi/ii|ei/i) + XnHmx/l - \c\ 2 \e n f n ). (33) 

Note that Eq. I|33|l expresses L\i/j) already in its Schmidt 
decomposition. 

The way to proceed now depends on the entanglement 
measure we choose. 



1. The renormalized measure of entanglement Ei 

In this subsection we will use the entanglement mea- 
sure Ei(k\iP)) := E(tp)). Let us look at the state \ip(ci)} 
of i|32|) with ci = 1/V2- This is obviously a maximally 
entangled state with Schmidt rank 2: Ei(\tp(ci))) = 
E(\ip(ci)}) = 1. The state lj3l?|) has this entanglement 
value only if Ai/ii = A„/i m , which fixes A n = Ai and 
Mm = Mi- (We remember that the /i's and z/s are or- 
dered.) So Da = Ai 1a and Db = Mi ^-b, which then 
leaves us with 



La = XiVaUa, 



L b =ixiV b Ub, (34) 



which are multiples of products of unitary matrices. This 
concludes the proof of Result for the renormalized 
entanglement measure. 



2. The probabilistic measure of entanglement E2 

We will now study the entanglement measure 
E2(c\tp}) = |c| 2 £'(|V')). Let us start with the physically 
most reasonable case, namely Ai/ii < 1. Assume now 
that we actually have A„/i m < 1. Then the entanglement 
of the maximally entangled state |-0(-^)) = -i=(|ai&i) + 

\a>nbm)) cannot be preserved, since this state would be 
mapped to the multiple of a state that is at most as en- 
tangled as the input state, but the length (or the proba- 



2 



< 1, 



bility of getting this state) would be 
leading to an even smaller actual value of the probabilis- 
tic measure of entanglement. Therefore we would have 
in this case Ai = A n = 1 = Mi = Mm, and L is unitary. 

Since L is of full rank, it is invertible and L~ l is also 
entanglement preserving, so that similar arguments can 
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be used to deal with the case X n jj, m > 1. The only re- 
maining case is then Ai/ii > 1 > A„/i m , and we deal 
with that in Appendix^! since it is rather technical, and 
moreover this case involves the rather unreasonable ar- 
gument of the increase in norm for some states. 



III. CONCLUSION 

Summarizing, we have investigated the properties of 
maps that preserve either entropy in single quantum sys- 
tems or entanglement in composite quantum systems. 
The first part of the paper demonstrates that the evo- 
lutions that respect the second law of thermodynamics 
and are probabilistically linear, are just the ones which 
are postulated in Quantum Mechanics: linear (on super- 
positions of vectors) and unitary. In the second part of 
the paper we have shown, that if from a linear map L (lin- 
ear on superpositions of vectors), we additionally demand 
preservation of entanglement in a very reasonable form 
(precisely that product states are mapped to product 
states and entangled states to entangled ones), then one 
immediately obtains preservation of the Schmidt rank. 
From this fact we have concluded that every reasonable 
entanglement preserving map has to be either local or 
local after application of a swap operation. In a further 
step, we have studied the consequences that two special 
extensions of the unique asymptotic entanglement mea- 
sure for pure bipartite states, induce on the two local 
parts of L. We were able to show that for the proba- 
bilistic measure of entanglement Ei , the linear map L is 
in fact a local unitary, or a local unitary after the appli- 
cation of the swap. And in the case of the renormalized 
measure of entanglement E\ , the map is either a multiple 
of a local unitary, or a multiple of a local unitary after 
the application of the swap. But, in this case (of E±), we 
can drop the factor, since we renormalize the outcome, 
for the evaluation of the entanglement. Furthermore, we 
have shown along the way, that the addition of the local 
ancillas does not increase the set of allowed maps since 
the space mapped into has to remain isomorphic to the 
composite Hilbert space H.ab- 

It is tempting to think of our results in the following 
way: entropy is, of course, preserved in unitary dynam- 
ics. However, for systems whose dynamics is not assumed 
to be quantum, entropy preservation can be seen as the 
second law of thermodynamics (for reversible processes) . 
On the other hand, for composite systems, it is clear 
that entanglement is preserved by local unitary dynam- 
ics and (complete) exchange of systems. However, for 
composite systems on which we do not put any dynam- 
ical assumption, preservation of entanglement has been 
proposed as an equivalent second law of thermodynam- 
ics for composite systems (see e.g. |38|). The premises 
we use could, therefore, be regarded as second laws for 
single and composite systems plus some forms of linearity 
on the evolution. And under these assumptions we have 
shown that the evolution must be unitary. 
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APPENDIX A 

In this Appendix, we argue that the assumption of 
preservation of disorder in reversible processes (item(i) of 
II A|) can coexist with statistical mechanics, in particular 
with the Boltzmann H theorem. 

Suppose that there are N classical particles of a gas 
in a box, such that the probability density at position 
r, momentum p, and time t, is given by the distribution 
function f(r,p,t). Then if the particles follow classical 
mechanical equations, the quantity 

Hit) = - J drdpf(r, p, t) log 2 f(r, p, t) 

can be shown to be an increasing function (under cer- 
tain assumptions). This is the Boltzmann H theorem 
[32ll33j| . Identifying H(t) with the disorder of the system, 
it may seem that the Boltzmann H theorem is against 
the hypothesis of preservation of disorder. However the 
H theorem is true in the limit of "molecular chaos" [34| , 
which means that in a volume element dr around the 
position r, the probability of finding two particles with 
momenta p\ and pi are independent, at any instant of 
time t, so that the probability of finding them simultane- 
ously is the product f(r,px,t)f(r,p2,t) of the individual 
probabilities. It is known that under such an assumption 
of factorization, an otherwise closed system turns into an 
open system, leading to irreversibility in the process (see 
e.g. H|). 

Our considerations are however for reversible pro- 
cesses, and so the H theorem is not relevant for our 
purposes. Remaining still within the realm of classical 
mechanics, the Liouville density D of a system of N (clas- 
sical) particles in the 6iV-dimensional phase space can be 
used to define the disorder of the N particle gas as the 
Gibbs entropy 

H G (t)=- J dTDlog 2 D, 

where dT is an element of the phase space. Using Liou- 
ville's theorem (stating that D is a constant of motion) 
|3t| , it is straightforward to show that the Gibbs entropy 
is a constant of (classical) motion. 

Similar considerations hold in the quantum domain. 
It is known |ll| (see also [37J) that a quantum system 
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approaches its equilibrium state through multiple colli- 
sions with other systems in equilibrium. In the process of 
reaching this equilibrium, the von Neumann entropy of 
the system increases, compatible with the average energy 
of the system. In the absence of any such interactions 
with other external systems, the von Neumann entropy 
of the system remains constant, as a consequence of the 
unitarity of its evolution. For an interacting system, the 
transformation of the system along with the heat bath 
is again unitary, and hence the entropy of the system 
plus the heat bath is again preserved. However, when 
we consider the system only, we see an increase of en- 
tropy. Such a system is open, and leads, quite generally, 
to irreversibility. 



APPENDIX B 

In this appendix we close the gap in the argument con- 
cerning the probabilistic measure of entanglement. We 
have to deal with the case when for the Xi/ii and \ n \x m 
of Eq. H8ri|) . the inequality \\\x\ > 1 > A„/i m holds. Due 
to continuity of the lengths of the vectors, there exists 
a C2 such that the the normalized vector \ty{caf) (of the 
form in Eq. is mapped by L to a normalized vector 

(in the form in Eq. (|33[l ). i.e. 



\tp(c 2 )) := c 2 |ai&x) + cj\a n b m ) 

±>l\x iyi ) + y/l-l 2 \x n y m ). (Bl) 

Due to preservation of entanglement by L the Schmidt 
values of the two vectors have to be preserved, and we 
have that either c^Af/if = c| and (1 — c|)A n /x m = (1 — 



c|), which immediately results in Ai/ii = 1 = A n /i m , or 
c|A?/ii = 1 — c| and (1 - c$)\ n fi ri 
equations can be combined to 



c|. The last two 



Ai/LiiA n /i m — 1. 



(B2) 



Now, since L is of full rank, every vector is in the range 
of L. So there exists a C3 such that 



\i>(c 3 )) := c 3 |ai&i) + <Jl- c\\a n b m ) 



1 



— ► \/ E ( c 3)-^(\xiyi) 



(B3) 



where E(c 3 ) is the entanglement of (^(03)), which only 
depends on C3. We therefore have the requirement 



Aipi = 



24 • 



E(c 3 ) 



2(1 



4Y 



(B4) 



Using the fact that the product of the two terms has to 
equal 1, as shown before, we have that 



E(c 3 ) 



= 2, 



(B5) 



which can easily be seen to have the only positive solution 
C3 = -j=. But this corresponds to the case where it is a 
maximally entangled state that is taken to a maximally 
entangled state, and Ai/ii = 1 as well as A„/i m = 1. So 
one can take Xi = 1 and fij — 1, which shows by Eq. 
(|34|l that La and Lb are actually unitary maps. So L is 
the product of two local unitaries, and not the multiple 
thereof, like in the case of the renormalized measure of 
entanglement. 
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